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General properties of the effective conductivity cr e of planar isotropic randomly 
inhomogeneous two-phase self-dual systems are investigated. A new approach 
for finding out cr e of random systems based on a duality, a series expansion in 
the inhomogeneous parameter z and additional assumptions, is proposed. Two 
new approximate expressions for cr e at arbitrary values of phase concentrations 
are found. They satisfy all necessary inequalities, symmetries, including a dual 
one, and reproduce known results in various limiting cases. Two corresponding 
models with different inhomogeneity structures, whose <r e coincide with these 
expressions, are constructed. First model describes systems with a finite maxi- 
mal characteristic scale of the inhomogeneities. In this model <r e is a solution 
of the approximate functional equation, generalizing the duality relation. The 
second model is constructed from squares with random layered structure. The 
difference of <r e for these models means a nonuniversality of the effective con- 
ductivity even for binary random self-dual systems. The first explicit expression 
for cr e can be used also for approximate description of various inhomogeneous 
systems with compact inclusions of the second phase. The percolation problem 
of these models is briefly discussed. 



PACS: 72.80.Ng, 72.80.Tm, 73.61.-r 

1. Introduction 

The electrical transport properties of the classical inhomogeneous condensed 
media, being a particular case of the transport properties in disordered systems, 
have a great practical interest. For this reason they have been intensively studied 
theoretically as well as experimentally (see, for example, 012)- An importance 
of this knowledge has revived last time. It was established that many real 
systems, which have been considered to be homogeneous, can have on small 
scales (from microscopic through mesoscopic till macroscopic) inhomogeneous 
or even heterophase structures. In particular, such situation takes place in 
the oxide materials, which are now actively studied |3J- Analogous situation 
takes place in many real examples of the disordered systems (for example @]). 
For this reason the observed effective conductivity a e of these systems can 
differ from that of homogeneous media. Thus one needs a theory adequately 
describing conducting properties of such inhomogeneous media. Unfortunately, 
the theory of strongly inhomogeneous systems is not well developed due to a 
difficulty of explicit averaging over strong inhomogeneous fluctuations. The 
various perturbation, approximate and variation theories have been developed 
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and many interesting results have been obtained so far (see, for example, jS] 
and a recent review [§]). 

In the conducting properties of the inhomogeneous media there is one old 
problem about the effective conductivity a e of inhomogeneous (regular, irreg- 
ular or random) heterophase systems, containing N(N > 2) various phases 
with different partial conductivities <Ji (i — 1,2, ...,N). In general, the inho- 
mogeneities can be very different, changing from regular through weakly and 
moderately irregular till strongly irregular ones. Strong inhomogeneous fluctu- 
ations can be connected with large differences of the partial conductivities as 
well as with strong irregularities of these phases. The latter situation is very 
similar to the case of random inhomogeneous medium and in some cases can be 
identified with it. 

In principle, the regular and weakly irregular cases can be considered sys- 
tematically in framework of classical potential theory |Jj , while an investigation 
of strong and random inhomogeneous systems needs some approaches connected 
with the averaging over the random regions with different conductivities |S] . 

The basic equation, describing a stationary electric current in random sys- 
tems, is the Ohm law, connecting the local current j(r) and the local electric 
field E(r), 

j(r) = <r(r) E(r), (1) 

where cr(r) is a local conductivity. We will assume that the carriers in all phases 
have the charges of the same sign. Then one can suppose that a(r) > 0. The 
equation (1) must be supplemented with the corresponding boundary conditions 
on the boundaries of two phases 

jl=jl Ej=El (2) 

where n,t denote the normal and tangent components and 1,2 correspond to 
different phases. The vector fields E(r) and j(r) satisfy the following equations 
inside each phase 

VxE(r)=0, V-j(r) = 0. (3) 

They mean that E(r) and j(r) are potential and divergenceless fields respec- 
tively. 

The effective conductivity a e of the isotropic randomly and/or strongly 
inhomogeneous medium can be defined as a proportionality coefficient between 
mean values of j and E averaged over the area of the system 

J = TeE, (4) 

}=fi(v)dr/S, E = jE(r)dr/S (5) 

where S is an area of the system. 

Since it is very difficult to make an averaging of (4,5) over an area of the 
system explicitly for strongly irregular or random systems it is usually supposed 
that for systems with such c(r) (or with a fixed random realization of cr(r) the 
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averaging of (4,5) over the area of the system can be replaced by the averaging 
over some ensemble of random functions er(r) with a functional weight V{cr}, 
which may, in general, depend also on temperature T 

J drF{a(r))/S = J Va{r)F{a(r))V{a} = < F{a(r)) > . (6) 

In this paper we confine ourselves with an investigation of the a e dependence 
on <7j and assume that T is fixed. 

According to the general principles, a behaviour of the effective conductivity 
of random systems must depend on statistical properties of the random conduc- 
tivity cr(r) encoded in its correlations. In case of random heterophase systems 
one usually supposes that cr(r) at different space points are statistically inde- 
pendent and its moments in each space point r are defined by one distribution 
function of partial conductivities -P(ct), which does not depend on r, 
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P(cr) = ^x;<5(ct - o-j), y^ Xj = l, / P(a)dcr = 1, 
i=l l ■* 



(7) 



where Xi is the concentration of the i-th phase. 

This means that a mean conductivity <C cr 3> (and all its higher moments) 
can be represented in two equivalent forms 

<cr>=y cr(r)dr/S = J crP{a)dcr = (cr). (8) 

where (. . .) denotes the average with a distribution function P(cr) Then the 

bare phenomenological correlation function of conductivities Co(r) takes the 
simple form 

C (r, r') =< a{r)a(r') >= (a) 2 + (a 2 )S(r - r') (9) 

As follows from (9) Co(r, r') in this form does not already describe a detailed 
spatial structure of random inhomogeneities, but contains only their integral 
characteristics- the concentrations of phases. In this case a e can be represented 
as a function of partial conductivities <Ti and concentrations xi only 

cr e = a e {ai,xi\a 2l x 2 \:-\cr N ,XN) (10) 

Due to the linearity of the defining equations this function must be a homoge- 
neous function of degree one of <Tj , i = 1, N. It has also satisfy the inequalities 
following from general principles (see, for example, 
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<a e <(a). (11) 



For practical purposes and for a description of experimental data one needs usu- 
ally some explicit and compact (though, maybe, approximate) formulas for effec- 
tive conductivity applicable in a wide range of parameters. A general expression 
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for a e of isotropic media at arbitrary concentrations Xi can be obtained only 
in case of weakly inhomogeneous media, when the conductivity fluctuations Aa 
( Aa = a — (a) , are much smaller than a mean conductivity (a) pQ 



where D is a dimension of the system. This formula expresses a e only through 
two first moments of random variable a. The further theoretical results in 
this problem for random inhomogeneous systems are related with the effective 
medium approximation (EMA) [5] and the random resistor networks (RRN) 
models of these systems [2], since they allow to fulfil an averaging procedure 
explicitly. The RRN models approximate the continuous inhomogeneous sys- 
tems (IS) and can be obtained as a result of some their discretization (though, 
in general, this procedure is not exact and unambiguous). It was shown for the 
RRN models firstly by numerical simulations 2 and later by series expansion 
on moments of 5a jlOllllj that the EMA, dealing ad hoc with weakly IS, works 
well even for moderately inhomogeneous RRN models. 

It is widely believed that the EMA, due to its good description of RRN 
type models, is also a good approximation for continuous random IS. As is 
known the EMA for different isotropic RRN models (for example, on isotropic 
regular lattices) gives different expressions for a e 0- Thus a e in the EMA is 
not an universal function even for isotropic RRN models. At the same time, 
not every continuous random IS can be described by the discrete RRN models, 
since their reduction to the RRN type models can change some properties. 
Moreover, the continuous random IS can have more complicate local structures 
(see below) . Then a e of such continuous systems can differ in higher order terms 
in 5a (starting with the 3-rd order, because the two first terms are universal 
and are determined by the formula (12)) from the EMA. For example, it is 
not obvious that any system, containing random "resistors" and "conductors" 
(i.e. the elements inversed to each other) simultaneously, must have the same 
behaviour as a "pure" RRN model. For a description of continuous IS it is 
very desirable to have any method, basing on the first principles and general 
properties, which is not connected with the RRN type models and/or with 
the EM approximation. Fortunately, for two-dimensional systems, due to the 
existence of the exact duality relation ^1 E|, one can attempt, basing on 
this relation, to find other possible approximations different from the EMA 
and applicable for systems with other inhomogeneity structures. In this paper 
we will show for the binary systems ( N = 2 ) how starting with a duality 
relation (DR) and using some additional assumptions one can find two explicit 
approximate expressions for a e , differing from the EMA. We will also construct 
two corresponding models, having their effective conductivities just of these two 
forms, and discuss their properties and possible applications. A case of arbitrary 
number of phases N will be presented in the subsequent paper. 




(12) 



2. Reciprocity and duality relations and symmetries 
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We begin with a consideration of main general properties of a e , because 
in the existing literature there are some misunderstandings about them. For 
systems with 5 -correlated inhomogeneities <r e can be reduced to a function, 
depending on the pairs of partial parameters: conductivities er,: and phase 
concentrations Xi (0 < Xi < 1, i — 1, 2, J^i x i = 1) 

cr e = (j e (a 1 ,x 1 \a 2 ,x 2 ). (13) 

Note, that in general case the order of the arguments pairs is important, since 
the phases can have different geometrical and/or statistical properties. In two 
dimensions there is a reciprocity between divergenceless and curl free vector 
fields v and v 

Vi = e lk v k , diVi = e ik diV k = 0. (14) 

This relation is local and takes place in each point. It means that the reciprocal 
vector field is the initial field turned on it/2. For this reason for continuous 
inhomogeneous (and isotropic) media the reciprocal vector field feels the same 
inhomogeneities of the initial system. Inside any phase one can connect with 
the reciprocal vector fields j and E a reciprocal conducting system, such that 
j r ~ E and E r ~ j, and for which the corresponding Ohm law takes place 

j r = <rE. r , a = al/a, (15) 

where <To is some arbitrary constant with dimension of conductivity. It means 
that a is inversed to a and the coordinate axis of the reciprocal system are 
turned on tt/2. Note that the corresponding phase boundary conditions (2) 
are also fulfilled, since the normal vectors transform into tangent ones and vice 
versa. The effective conductivities of the initial medium and its reciprocal must 
satisfy after averaging the same functional relation as in (15), which can be 
called the reciprocity relation (RR) and is usually written in the following form 

Ve(vi,Xl\<T2,X 2 )Ve{vi/vi,X 1 \<jl/(T2,X2) = (16) 

In this general form it have been obtained firstly by Keller ^2] an d indepen- 
dently by Dykhne (T^|. The general RR (16) takes place for all two-dimensional 
binary systems: for RRN models and for regular and irregular media, as well 
as for random ones. But the general RR is not so useful, since it connects 
the effective conductivities of different systems. In order to obtain more in- 
formation about <r e one needs to know how these effective conductivities are 
related. The most interesting case is when a reciprocal system has the same 
statistical or/and geometrical properties. In this case the effective conductiv- 
ity of a reciprocal system must coincide with the effective conductivity of the 
initial one. Such systems can be called the self-reciprocal ones (we will use a 
name self-dual for more symmetrical systems, see below). One can see from 
the above description of the reciprocal medium that for isotropic continuous 
random inhomogeneous medium the reciprocal medium will have the same type 
of inhomogeneities and after averaging its effective conductivity a e must have 
the same functional form as the effective conductivity of the initial medium, but 



5 



with inversed partial conductivities. Thus they have satisfy the self-reciprocity 
relation (SRR) 

a e (a 1 ,xi\a 2l x 2 )a e (al/(Ji,xi\al/a 2l x 2 ) = <J%. (17) 

At the same time it is worth to note here that in case of discrete RRN type 
models the reciprocal network can differ from the initial one and consequently a 
functional form of the effective conductivity of the reciprocal network will differ 
from that of initial network. In particular, the effective conductivities of RRN 
models on isotropic regular lattices (let us call them random resistor lattice 
(RRL) models), triangle and honeycomb, turn out different and must satisfy, 
due to a mutual reciprocity of these lattices, only the R-relation 

a t e r (a 1 ,x 1 \a 2 ,x 2 )cr^ on {(Tl/a l7 x 1 \al/(T2,X2) = (18) 

This is a payment for an ambiguity of the discretization procedure of the con- 
tinuous medium. 

The self-reciprocal systems include the continuous media and RRN model on 
self-reciprocal networks. For self-reciprocal systems the SRR can be rewritten 
also in the form 

a e (ai 1 xi\a 2 ,x 2 ) = al/a e {al/(T ll xi\(jl/(j 2 ,x 2 ). (19) 

This form has a simple physical sense, connected with two possibilities of a 
writing the Ohm law in terms of conductivities a (as (15)) and in terms of 
resistivities p = 1/cr : an effective conductivity of a self-reciprocal system as a 
function of partial conductivities must coincide with an inversed resistivity as 
a function of corresponding partial resistivities (i.e. inversed partial conductiv- 
ities). For this reason the relation (19) can be also called the inversion relation 
(IR). 

Though now the SRR (17) restricts a possible functional form of a e , in 
a general case it does not yet permit to find the exact values for a e even at 
equal phase concentrations x\ = x 2 = 1/2 and (Tq = <j\<j 2 (the latter choice 
turns the inversion a — > cxq /a into permutation <ji <-> 02 ). To be able to find 
exact value of a e one needs additional symmetry property of systems connected 
with the geometrical and/or statistical properties of their partial phases. To 
see various possibilities, appearing for inhomogeneous media, let us consider a 
simple, but very instructive, case of regular inhomogeneous system consisting 
from one (white) phase with regular round inclusions of the second (black) 
phase, the Rayleigh model. This system is self-reciprocal in the above sense 
and its effective conductivity satisfies the SRR (and IR) (17) (52)- However, 
the geometrical properties of these phases are different. The four pictures of 
this model at the interchanged phases and at the adjoint concentrations x and 
1 — x of black phase are represented on Fig.l. As follows from Fig.l, all these 
examples of the model are different and, consequently, the corresponding a e 
are also different. Thus, though the SRR takes place for this model, it does not 
give any equation for a e . Such equation can be obtained only when a e is a 
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(a) 



(b) 



(c) (d) 
Figure 1: (a) A part of the model with regularly situated round inclusions at the 
concentration x , (b) a part of the same model at adjoint concentration 1 — x , (c) 
the same model at concentration x but with the interchanged phases, (d) the same 
model with interchanged phases at adjoint concentration 1 — x . 



symmetrical under simultaneous interchange of both parameters of phases 

(T e {a\,xi\cr-2,X2) = a e (a2,X2\ai,xi). (20) 

This property can take place as for random systems as for some regular systems. 
The typical examples of such systems are the random mixtures and RRN models 
with the probability distribution function P(a) (7), and a regular "chess-board" 
type system. One can see from Fig.l that the self-reciprocal Rayleigh model 
does not have this symmetry, but this symmetry appears as an approximate 
symmetry, when a concentration of one phase is small (then fig. la (fig. lb) is 
similar to fig. Id (fig.lc)). 

We will call symmetrical self-reciprocal systems as self-dual (SD) systems. 
The corresponding SRR takes at CTq = u\a 2 the following form 

cr e (ci, xi\a 2 , x 2 )a e (ai, x 2 \a 2 , X\) = a\a 2 . (21) 

For SD systems one can obtain exact equation for a e . It takes place only at 
equal phase concentrations Xi = 1/2 and at Uq = aia 2 . and gives the well- 
known exact Keller-Dykhne value |12l IT5| 

o- e (oi, l/2\a 2 , 1/2) = (22) 
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This remarkable formula is very simple and universal in a sense that it takes 
place for all two-dimensional inhomogeneous self-dual systems as for regular as 
for random. 

For binary systems one can write a e in a more short form as a function of 
its independent parameters only (here x is a concentration of the first phase) 

a e = a e (x,ai,a 2 ). (23) 

Then, for isotropic SD systems with unequal phase concentrations (21) can 
be written in the form, which connects the effective conductivities at adjoint 
concentrations x and 1 — x or in terms of a new variable e = x— 1/2 (— 1/2 < 
6 < 1/2 ) at e and — e 

a e (x,ai,a 2 )<J e (l - x,ai,a 2 ) = <7\<J2 (24) 

Eq.(24) means that a product of the effective conductivities at adjoint concen- 
trations is an invariant. Due to this relation one can consider <j e only in the 
regions x > 1/2 ( e > ) or x < 1/2 ( e < ). Just the relations (21) and (24) 
are usually named as the duality relation (DR). Below we will consider mainly 
systems, satisfying the DR (21) or (24) as the exact functional relation. It fol- 
lows also from (24) that (for > er 2 ) than faster is a e increasing at x > 1/2, 
than slower it increase at x < 1/2. 

The binary case assumes further simplifications. It is convenient to use 
instead of <7j,i = 1,2 another combinations of Cj : a± — (<ti ± <72)/2. Then, 
the effective conductivity can be represented in the form, symmetrical relatively 
to both phases, 

CT e (e, er+, cr-) = cr+/(e, cr-/er+) = <r+/(e, z), (25) 

where a new dimensionless variable z = er_ / 'cr + , (— 1 < .z < 1), characterizing 
an inhomogeneity of the system, is introduced. This choice of variables is con- 
venient for a study of a e in the limit, when one partial conductivity becomes 
small or goes to 0. 

The effective conductivity <r e (e, a +1 (7_) and a function f(e,z) must have 
the next boundary values 

cr e (l/2,er + ,(7_) = CTi, C7 e (-l/2, cr+, <7_) = <T 2 , 

f{l/2,z) = l + z, f(-l/2,z) = l-z, /(e,0) = l. (25') 
The duality relation takes in these variables the form 

f(e,z)f(-e,z) = l-z\ (26) 
from which one obtains at critical concentration e = the exact value 

/(O, z) = Vl-z 2 - (26') 
The self-reciprocity relation in these terms has the form 

f(e,z)f(e,-z) = l-z\ (27) 
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and cannot give a nontrivial exact value, since /(e,0) = 1. The symmetry (20) 
in the new variables means that 



f(e,z) = f(-e,-z), f(-e,z) = f(e,-z). (28) 

It follows from (28) that the even ( f s ) and odd ( f a ) parts of /(e, z) relatively 
to e coincide with the even ( f s ) and odd ( f a ) parts of f(e, z) relatively to 
z. The simplest way to satisfy these properties corresponds to the following 
functional form of f{e,z) 

/(e,z) = /(ez,e 2 ,z 2 ) (29) 

One can conclude from (29) that in this case an expansion of f(e, z) near the 
point e = z — does not contain terms linear in e and z separately. The odd 
part f a can be represented in the form 

/ Q (e,z) = 2ez$(e,z) (30) 

where $ is an even function of e and z (the coefficient 2 in front of ez is 
chosen for further convenience). 

The DR (24) in general case is not enough for a complete determination of 
/ , but it poses some constraint on the possible form of /. There are various 
ways of its choice. The most known example, using the partition of / into 
symmetric and antisymmetric parts / = / s + f a , connects f s and f a 

/ 2 -f a =l- z\ (31) 

Since this form of the DR is associated with "a parity" partition of function / 
one can call it "a parity representation" of the DR. It means that f a and /„ 
considered at fixed z as the functions of e satisfy to hyperbolic relation with a 
constant depending on z. The relation (31) allows to express f(e, z) through 
its even or odd parts 

f(e, z) = f a + V/ 2 + 1 - z 2 = f s ± y/f? -l + z 2 - (32) 

For this reason it is enough to know only one of these two parts. Usually 
one prefers to choose an antisymmetric part as more simple. Using the first 
approximation in z for f a in the weakly inhomogeneous case, which follows 
from the mean conductivity 

(a)=a+(l + 2ez), f a (e,z) = 2ez. (33) 

and substituting (33) into (32), one obtains the effective conductivity 



a e (e,z) = a + 2ez + y/(2ez) 2 + 1 - z 2 , (34) 

which coincides with a e of self-dual square RRL model in the EMA [2]. We 
will refer the expression (34) as the EMA expression for effective conductivity of 
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self-dual systems. For self-reciprocal systems there are also relations analogous 
to (31,32) with a change f s , f a — > f s , /°, but in this case the relations (29,30)do 
not take place. For this reason the EMA cannot be a good approximation for 
SR systems. 

For a checking the symmetry properties and a comparison of different expres- 
sions for <j e it is useful to consider also a e of RRL models on other isotropic 
lattices, triangle and honeycomb. The EMA expression for a e of the binary 
lattice RRN models is [2] 

ae = 7^2 ( D + v/[] 2 + 2(9-2)^^} , 

0= [(fx- l)*i + (f(l-*) -l)o- 2 ] (35) 

where q is a number of bonds coming in (or outgoing from) one lattice site. In 
terms of z and e it has the following form 



q-A 



,1/2 



q-2 { V 2 1 

[] = [(q/2) 2 + q(q - 4)ez + (q 2 e 2 - 2(q - 2))z 2 ] (36) 

which has the necessary symmetry property (29). One can check that the cor- 
responding ai tr ^ (q=6) and ai hon ^ (q=3) satisfy also the RR 

As it was remarked earlier, some properties of RRN models differ from the 
properties of continuous media. For this reason, one can hope that the other 
possible representations for / compatible with the DR can produce other ap- 
proximate expressions for er e , which will turn out more suitable for models with 
other inhomogeneity structures. Another reason to analyse the duality relation 
is connected with the known fact, that usually the systems with such symmetry 
have some additional, hidden, properties, permitting to obtain more informa- 
tion about function under consideration. Moreover, in some simple cases they 
can help to solve problem exactly (see, for example, |15p. Bearing this in mind, 
we will investigate various consequences of the duality relation in more details. 



3. Duality, symmetries and perturbation expansions 

We start with the series expansion of /(e, z) on z. For every fixed z (0 < 
z < 1) a function / must be a monotonous function of e. Since a homogeneous 
limit z — ► is a regular point of / it will be very useful to consider a series 
expansion of / in z 

oo 

/(e,z)-^A.(e)z fe /fc!, (37) 
o 
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where due to boundary conditions (25') 

f = 1, /i(e) = 2e. (38) 

It is worth to note here that: 1) similar expansion was already considered up 
to third order (see, for example, Hi ) and 2) the expansion (38) is connected 
with the weak disorder expansion for RRN models in moments of fluctuations 
da We will show that our expansion is much simpler, since it treats the 

variables z and e independently. From this point of view the expansion on 
moments of 5a is an expansion on rather complicated functions of z and e. 
Substituting the expansion (37) into (26) one obtains the following results: 

(1) in the second order in z it reproduces the universal formula (12), thus 
the latter can be considered as a consequence of the duality relation; 

(2) in general case, there are the recurrent relations between f 2 k and f^k-Xi 
corresponding to the connection (31) and expressing the even terms through the 
odd ones; 

(3) /2fc+i(e) are odd polynomials of e of degree 2/c + l and /2fc(e) are even 
polynomials of e of degree 2k in agreement with (29). 

Taking into account boundary conditions (25'), the exact value (26') and the 
symmetry property (29), one can show that the coefficients fk must have the 
next form 

/ 2 fe+i(e) = e(l-4e 2 ) 52fc _ 2 ( e ), k > 1, 

f2k(e) = (l-4e 2 )/i 2fc _ 2 (e), k > 1, (39) 

where <?2fc-2(e) and 7i2fe-2(e) are some even polynomials of the corresponding 
degree, and free terms of /i2fe-2(e) coincide with the coefficients in the expansion 
of the exact value (26') 

\]\ - z 2 = 1 - z 2 /2 - z 4 /8 - z 6 /16 - z 8 /128 - z 10 /256 + ... (40) 

The parity representation gives the relations, expressing / 2 fc with even indices 
through fk with lower indices. For example, 

h = 4/i/ 3 - 3/1 = (1 - 4e 2 )[(8.g + 12)e 2 - 3], (41) 

h = 6/1/5 - 15/ 2 /4 + IO/3 2 . (42) 

It follows from (39) that /a is completely determined up to overall factor 
number go- Since fi is determined through lower fk (fc — 1,2,3) it is also 
determined by the coefficient go. In the EM approximation the series (21) is 
defined completely and takes very simple form, since in this case all g 2 fc-2 = 
0(fc>l) and / 2fe (e)~(l-4e 2 ) fe . 

It is useful also to write out the series expansions for other 2D isotropic RRL 
models in the EM approximation. The expansion of (36) in z up to the 6-th 
order terms gives 

L, 2 2 2 4e(l_4e 2 K8_6 i +a 3 
a e = a+ < 1 + 2ez (1 - 4e )z H Tl rr z a 

I q <r(<7-2) 
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4(l-4e 2 )(< ? -2 + 2e 2 (20-10g + g 2 )) 4 



' ' ' 4g2) ^ - 6g + g 2 )(3(g - 2) + 2e 2 (28 - Uq + q 2 ))z 5 - 



1,21 

\ 2 I rt 2 / i i n I iio n Q/i „2 | o „3 



16(1 5 4f2) ((«? - 2) 2 + 2e 2 (-112 + 112«? - 34 g 2 + 3g 



Q 

+2e 4 (672 - 672g + 224q 2 - 28g 3 + q 4 )) z 6 } (43) 

It has the necessary structure following from the general properties. Note that 
the two first terms coincide with the universal (12) only for self-dual square 
RRL. One can see also from (43) that the higher odd terms (starting from the 
3-rd order term) depend on the lattice structure and do not vanish except the 
self-dual case of the square lattice (where this takes place due to the factor 
(8-6<? + <7 2 )). 

In order to find higher order terms for square RRL model one may compare 
the expansion (37) with known results in scries expansion in powers of conduc- 
tivity fluctuations 8k = ((Sa) k ) 11 . Its coefficients are found up to 6-th order 
and it was shown that for 2D square RRL model, due to the duality relation, 
they begin to differ from the corresponding EMA coefficients only in the 5-th 
order 

8 2 , 8 3 84 : 8 5 + J18283 



(44) 



2 4 8 16 

[S 6 + (3Ji - 1)S 2 8 4 + (2Ji - 1)81 - (gJx - 2) 
32 



where J\ =~ 1.092958179 is some lattice integral of the square lattice Green 
function. The EM approximation corresponds to J% = 1. 

For a comparison of the general expansion (44) with (37) one needs the 
explicit expressions for the higher moments of 8a . After some algebra one can 
show that the moments of binary systems with the distribution function (7) 
have in terms of variables a— and x (or z and e ) the following form 

((Aa) k ) = x(l - x)P k - 2 {x){cj-) k - (a + ) fe (l - 4e 2 )/4P fe _ 2 (e)z fc , 

Pk-2(x) = {l-x) k - l + {-l) k x k -\ 

P fe _ 2 (e) = (1/2 - e)*- 1 + (-l) fe (e + l/^ 1 (45) 

One can see from (45) that z -dependence of the moments is factorised. More- 
over, it follows from the explicit expression for the polynomial Pu-2 that for 
odd k = 2m + 1 , due to the symmetry of this polynomial relative interchange 
x <-> 1 — x, it has additional general factor (1 — 2x) = —2e in accordance with 
(39). Thus, we see that the series expansion for square RRL model in 8u from 
has in terms of variables z and e the necessary structure, following from 
symmetry and duality. The obvious advantage of series expansion (38) is that it 



12 



treats the variables z and e separately. The explicit forms of the polynomials 
Pk-2 at k < 6 are 

(A 3 ) = 4e(l - 4e 2 )a 3 _ , (A 4 ) = 4(1 - 4e 2 ) (1 - 3x + 3x 2 )a 4 _ , 

(A 5 ) = 4e(l - 4e 2 )(l - 2x + 2x 2 )a 5 _, 

(A 6 ) = 2 4 (1 - 4e 2 )(l - hx + lOx 2 - 10x 3 + 5x 4 )a 6 _. (46) 

Using the equality (33) for (a) and expanding the terms of the series (44) in 
z, one can obtain the series (37) for the square RRL model up to 6-th order. 
In particular, in the EMA one obtains the corresponding series (43). Since the 
series for square RRL model begin to differ from the corresponding EMA in 
5-th order it is enough to find only 5-th and 6-th order terms. They are 

.g 2 = -30.7(1 -4e 2 ), h 4 = -45(1 - 4e 2 )[l - 4(1 - 2j)e 2 ], (47) 

where j = J\ — 1 = 0.092958179 <C 1. One can check that these formulas satisfy 
relations (41,42). 

Thus, we have obtained the series in z for <j e of self-dual square RRL 
model up to 6-th order. Since it coincide with the EMA up to 4-th order, 
one can say that they belong to the same class of functions with g = 0. As 
follows from (43) the effective conductivity for isotropic RRL models in the EM 
approximation depends on the lattice structure and can have go ^ 0. Then the 
next question appears: is the effective conductivity of the continuous self-dual 
systems universal and always belongs to the EMA class or there are SD systems 
with other functional forms of er e ? The answer on this question will be obtained 
in the subsequent sections. 

4. Exponential representation and new approximate expressions 

Thus we have seen that the arbitrariness of / of self-dual systems is strongly 
reduced by boundary conditions and the exact value (26'), and that the third 
and fourth order terms are determined only up to one constant g , which equals 
for the square RRL model and in the EM approximation. The question about 
the universality reduces in this particular case to the next one: is the equal- 
ity go = true for all self-dual random inhomogeneous systems or it depends 
on the structure of random inhomogeneities? If the value of g depends on 
random structure, then the second question appears: can any other additional 
information about the function / of these systems determine the constant g 
or even the whole function, as it was in case of the EMA? Since the EMA is 
connected with the parity representation (31) and a mean conductivity approx- 
imation (33) for odd part, one needs to know what kind of other representations 
for / can satisfy the duality relation (26). In order to answer on this question 
it is convenient in the case z ^ 1 to pass from / to / = //VT — I 2 . Then 

f(e, z)f(~e, z) = l = f(e, z)f(e, -z). (48) 
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The DR (48) assumes other possible functional forms of f(e,z). Taking into 
account (29), one can write out the exponential representation for / 

f(e,z)=exp(ez<j>(e,z)), (49) 

where (f>(e, z) is some even function of its arguments. It is easy to see that it 
automatically satisfies eq.(48). Now let us find out is it possible to define g 
from some additional information about function 4>. For this we consider two 
simple ansatzes for 0(e, z) , analogously to the parity representation (31) and 
the EM approximation. 

1) In the first case we will assume that <j>(e,z) depends only on z. This corre- 
sponds to the exponential dependence on concentration. Such dependence takes 
place often in disordered systems |17j . Expanding the corresponding functions 
/ in series one can check after some algebra that now all polynomial coefficients 
can be determined unambiguously! The series for f a up to 5-th order has the 
form 

2 z 2k+l 

fa = ^2f2k+i ^ - j-yy , /afc+i = e(l - 4e 2 )52fe-2, fc=l,2, 

00 = -1, 5 2 = -(ll + 4e 2 ). (50) 

Another way to see this is to apply the boundary conditions (25') directly to 
the function (49). Then one obtains the whole function 

0(z ) = l/ zm i±£ f (e)Z)= (i±£\\ (51) 

The series expansion of (51) reproduces exactly the corresponding expansion 
from (50). Note that now go ^ 0. The even part f s can be found from (41,42). 

2) In the second case we will assume that 0(e, z) depends only on the combi- 
nation ez. This corresponds to the dependence of er e only on mean conductivity 
and/or mean resistivity, since (c ±:L ) ~ (1 ± 2ez). Expanding the corresponding 
functions / in series one can see that again all polynomial coefficients can be 
determined completely. The coefficients for f a up to 5-th order are 

.90 = -3, .g 2 = -15(l + 12e 2 ) (52) 

In this case one can also apply the boundary conditions (25') directly to the 
function (49). Then one obtains 

1 l + 2ez ~ ( l + 2ez \ 1/2 
0(e ' z)= 2^ ln I-2^< ^ Z)= \T=2T Z ) ■ (53) 

The series expansion of (53) also coincides exactly with the corresponding ex- 
pansion represented in (52). 

It is worth to note that an analogous exponential representation can be 
obtained for self-reciprocal, but nonsymmetric systems. In this case only the 
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second equality in (48) is fulfilled and the exponential representation has the 
form 

f{e,z) = e3qp(z<j>(e,z)), (54) 

where <f>(e, z) is an odd function of z. The boundary conditions (25') cannot 
now determine it completely. They only give 

«4->-4 b (£f)- 

The simplest way to satisfy these boundary conditions is to assume that the 
interpolating function at arbitrary e has factorized dependencies on e and z 
(it is an analog of the assumption 1) in the SD case) 

^) = ^ln(l±£), (55) 

here a(e) is an arbitrary function except its boundary values 

a(±l/2) = ±1/2. 

For self-dual systems a(e) = e, and just this equality ensures the exact Kcllcr- 
Dykhne value (22,26') at e — 0. Formally, the SR relation admits the self-dual 
solution (52) in nonsymmetric case too, but, as we will show below in section 7, 
(52) works well only for e, belonging to small regions near the boundary values 
e = ±1/2 (or at small concentrations of one of the phases 1 — x <C 1 ). Its 

applicability at intermediate values of e depends on how close is a non self-dual 
system to the self-dual one. 

Thus we have shown that under some simple (but physically acceptable) 
additional assumptions about a function f(e, z) one can obtain the new ap- 
proximate (a mean field like) expressions for /. Their physical meaning and 
properties will be discussed in two subsequent sections where we will present 
two simple models having the effective conductivity just of two forms found 
above. Here we only want to note that it is very important for further progress 
in this direction to find any other explicit solutions. In general case, assuming 
an expansion of (j) in double series in z and e, 

°° z 2k e 21 
<Ke,z)= ^-fcTTT 

k,l=0 

one can show that now again the boundary conditions (25') cannot define all 
coefficients completely. For example, and fi contain one free parameter 
<Pw ■ 9o = 6(0io - !)■ 

5. Finite maximal scale averaging approximation 

In this section we show that the solution (52) can describe a e of inho- 
mogeneous systems with compact inclusions of one of the phases. Firstly, we 
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introduce a new composite method of a construction of random inhomogeneous 
systems. In framework of this approach one divides an averaging procedure on 
two steps, corresponding to the averaging on two different scales. Then, using 
this approximate averaging, which we will call as a finite maximal scale averag- 
ing approximation (FMSA), we will obtain an additional equation for effective 
conductivity of inhomogeneous systems with compact inclusions. This approxi- 
mation can be considered as some modification of the mean field approximation 
and it is applicable to the random inhomogeneous systems with a finite maxi- 
mal scale of the inhomogeneities. Below in this section we will use temporally, 
for clarity, the concentration variable x instead of e. Let us consider planar 
two-phase randomly inhomogeneous system. It is easy to see that, in general, 
its effective conductivity will depend on a scale I of a region over which an 
averaging is done. This takes place due to the possible existence of different 
characteristic scales in the inhomogeneous medium. In the most general case 
there will be a whole spectrum of these characteristic scales. This spectrum 
can change from discrete finite one till continuous infinite one and will define 
the inhomogeneity structure of the system. For this reason this spectrum can 
depend on concentration x. It is obvious that the effective conductivity of the 
system can depend on this spectrum as a whole as well as on which area an 
averaging is fulfilled over. Suppose, for the simplicity, that the randomly in- 
homogeneous structure of our system has a scale spectrum with the maximal 
scale l m {x), which is finite for all x in the region 1 > x > 1/2 (or 1 — x in 
the region < 1 — x < 1/2). Let us assume that we know an exact formula 
for a e (x,z) of this system, which is applicable from scales I > l m . It means 
that this formula for a e {x, z) takes place after the averaging over regions with 
a mean size I > l m and does not change for all larger scales I 3> l m . Now we 
deduce an approximate functional equation for this a e (x, z). To do this let us 
consider the auxiliary square lattice with the squares of length II ^> l m and 
suppose that they have the effective conductivities corresponding to different 
values of the concentrations x\ and X2 with equal probabilities p=l/2 (see 
Fig.2) 

After the averaging over the scales I > II one must obtain on much larger 
scales I ^> II the same effective conductivity, but corresponding to another con- 
centration x = (xi +X2)/2. This is possible due to the similar random structure 
of different squares and due to the conjectures that in this model: (1) there are 
only two maximal characteristic scales l m (xi) = h [i = 1,2) and the lattice 
square size II ^> h, l m (x) , (2) the averaging procedures over these scales do not 
correlate (or weakly correlate) between themselves. Thus, for a compatibility, 
all concentrations must be out of small region around critical concentration x c , 
where h or l m (x) can be very large. We will call further this set of the assump- 
tions the finite maximal scale averaging approximation (FMSA approximation) . 
It can be implemented for systems with compact inhomogeneous inclusions with 
finite l m . From the other side the effective conductivity of the auxiliary plane 
on scales I 3> II must be determined by the universal Keller - Dykhne formula 
(22). Thus we obtain the next functional equation for the effective conductivity, 
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(b) 



Figure 2: (a) An elementary square of the auxiliary lattice with compact inclusions, 
(b) the auxiliary lattice, the numbers 1,2 denotes squares with the corresponding 
concentrations x\ and X2 ■ 



connecting a e (x, z) at different concentrations, 

a e (x, z) = \fa e (xi, z)a e (x2, z), x = (x\ + x 2 )/2, (xi,x^x c ). (56) 

It must be supplemented by the boundary conditions (25'). This equation can 
be considered as a generalization of the duality relation (21), the latter being 
a particular case of (56) at x\ + x 2 = 1. It follows from (56) that, due to 
the exactness of the duality relation, it really works for all concentrations x, 
except, maybe, a small region near x — x c and z = 1 (this region corresponds 
to the singular region of the percolation problem, see also below a discussion of 
the percolation limit). One can easily check that a e in low concentration limit, 
being a linear function of x, satisfies (56). Moreover, one can find an exact 
solution of this equation. It has an exponential form with a linear function of 
x (or e ) 

<r e (x, z) = g\ cxp(aa; + 6), (57) 
where the constants a, b can be determined from the boundary conditions 

a = — 6, exp b = o 2 /&i ■ (57') 

Substituting these coefficients into (57) one obtains 

a e ( x , z )^a 1 (a 2 /a 1 f- x) =a + VT^( K \^) ■ (58) 

The solution (58) satisfies symmetry relation (20) and can be represented in the 
exponential form 

a e (x, z) = o+\J\ — z 2 exp ^eln | — , (58') 

which coincides exactly with the case 1) from the section 4 with <p(e, z) = 
i In . Its odd part has a form 

f a (e, z) = sinh(ez(/)(e, z)), (59) 
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and automatically satisfies (31). 

It is interesting to note that the form of the solution (58) means that in the 
considered approximation one has effectively an averaging of the log a since it 
can be represented as 

log a e = (log a) = x log a\ + (1 - x) log <j 2 • (60) 

This coincides with a note made in J3| for the case of equal concentrations 
x = 1/2 and with analogous results obtained later for disordered systems in the 
theory of weak localization ^H] . The solution (58) satisfies also the inequalities 
(11) due to the Holder inequalities 

(a- 1 )" 1 < a e (x, z) = of*?") < {<j). (61) 

One can check that (58) reproduces in the weakly inhomogeneous limit the 
universal Landau - Lifshitz expression (12). In the low concentration limit of 
the second phase it gives 

1 — z 

a e {x,z) = a x (1 + (1 - x) log- h ...), 1 - x < 1, (62) 

1 ~\~ z 

what coincides with the Maxwell-Garnett formula in the weakly inhomogeneous 
case pp. Note that the expansion (62) contains the coefficient logarithmically 
diverging in the limit \z\ — ► 1 (<T2 — > 0). Such behaviour of the coefficients 
denotes the existence of a singularity in this limit as it must be due to the 
percolation singularity (see below a discussion of this percolation limit). 

As follows from the arguments of the method used above, the FMSA can be 
applied also for approximate description not only of random systems, but also for 
any (regular as well as irregular) inhomogeneous models with compact inclusions 
of one phase into another. However, the explicit approximate expression (58) 
obtained for a e will work well in a wide range of parameters only for systems 
close to the self-dual ones (see the notes in the sections 4 and 7). For example, 
this expression describes a e of the Rayleigh model with strictly round inclusions 
not very well for intermediate concentrations (it gives smaller values) as it shows 
the comparison from section 7 (it can be expected also from Fig.l). But it will 
work much better when the shape of these inclusions is deformed so that the 
corresponding pictures of this model (Fig. la and Fig. Id, Fig. lb and Fig.lc) 
become similar and symmetry property (20) is approximately fulfilled. 

6. Random "parquet" model 

Now we will show that the existence of different functional forms for the ef- 
fective conductivity denotes a nonuniversal character of this value for self-dual 
two-phase systems with different random structures. We will construct another 
model of two-dimensional isotropic randomly inhomogeneous two-phase system, 
using the composition method introduced above, and find a mean field like ex- 
pression for its effective conductivity a e (x,z) in case of arbitrary phase con- 
centrations x. The derived formula satisfies again to all necessary symmetries, 
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(a) 



(b) 



Figure 3: (a) An elementary square of the model with a vertical orientation, the 
dotted regions denote layers of the second phase; (b) a lattice of the model, the small 
lines on the squares denote their orientations. 

including a self-dual one, and coincides with the example case 2) from section 
4, realizing the second variant, when f(x,z) depends (in this approximation) 
only on one combination of variables ez. This dependence is described by the 
function analytical at small values of this variable. 

Let us consider the following planar model. There is a simple square lattice 
with the squares consisting of a random layered (or striped) mixture of two 
isotropic conducting phases with constant conductivities <Ji,i — 1,2 and the 
corresponding concentrations x and 1 — x. A schematic picture of such square 
is given in Fig. 3. 

The layered structure of the squares means that the squares have some pre- 
ferred direction, for example along the layers. Let us suppose that the directions 
of different squares are randomly oriented (parallely or perpendicularly) relative 
to the external electric field, which is directed along x axis. In order for system 
to be isotropic the probabilities of the parallel and perpendicular orientations of 
squares must be equal or (what is the same) the concentrations of the squares 
with different orientations must be equal p|| = p± = 1/2. 

Such lattice can model a random system consisting from mixed phase re- 
gions, which can be roughly represented on the small macroscopic scales as 
randomly distributed plots with the effective "parallel" and "serial" connec- 
tions of the layered two-phase mixture (Fig. 3). The lines on the squares denote 
their orientations. 

This layered structure can appear, for example, on the intermediate scales 
when a random medium is formed as a result of the stirring of the two-phase 
mixture or as a result of appearance of randomly distributed regions of " striped 
phases" in oxide materials [HJ ■ The corresponding averaged parallel and perpen- 
dicular conductivities of squares m \ (x) and a j_ (x) are defined by well known 
formulas 

CT||(x) = (a) = xgi + (1 — x)gi = <T + (1 + 2ez), 
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Thus we have constructed the hierarchical random medium (in this case a two- 
level one). On the first level it consists from some regions (squares of two 
different types) of the random mixture of the two layered conducting phases 
with different conductivities o\ and 01 and arbitrary concentration. On the 
second level this medium is represented as a random parquet constructed from 
squares of two such types with different conductivities ai | and a± , depending 
nontrivially on concentration of the initial conducting phases, and randomly 
distributed with the same probabilities p t — 1/2 (Fig. 3). This representation 
allows us to divide the averaging process approximately into two steps (firstly 
averaging over each square and then averaging over the lattice of squares) and 
implement on the second step the exact formula (22). This can be considered as 
some modification of the FMSA approximation. As a result one obtains for the 
effective conductivity of the introduced random "parquet" model the following 
formula, which is applicable for arbitrary concentration 



a e (e,z) = V^X*- 1 )- 1 = (T+y/l - z 2 f(e, z), f(e,z) 



1 + 2ez 



1 - 2ez 



1/2 



(64) 



This function has all necessary properties and satisfies SD equation (26), in- 
equalities (11) and the symmetry (20). It coincides with the solution 2) from 
section 4 and has another possible functional form, automatically satisfying the 
duality relation (48) 

f(e,z) = B(e,z)/B(-e,z) (65) 

with a function B(e,z) = [l+2ez]', which depends only on the combination ez. 
Strictly speaking, this form is not independent, since (65) can be also represented 
in the exponential form 

f(e,z) = exp(<f>(e, z)) 
where (p(e, z) is an odd function of e 

cj>(e,z)=]nB(e,z)/B(-e,z). 

It is interesting to compare this formula with the known asymptotic formulas. 
Let us consider firstly its behaviour for small phase concentrations. 

(a) In case of small concentration of the first phase one gets 

/ 2xz \ 

<J e (x,z) ~ CT 2 ( 1 + — —5 j • ( 66 ) 

It follows from (66) that an addition of small part of the first higher conducting 
phase increases an effective conductivity of the system as it should be. 

(b) In the opposite case of small concentration of the second phase 1-iCl 
one obtains 

i.e. an addition of the phase with smaller conductivity decreases <J e ff- It is 
worth to note that both these expressions for arbitrary values of the conduc- 
tivities o i and (72 differ from the Maxwell-Garnett formula and coincide with 
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it only in the weakly inhomogeneous case i < 1. It must be not surprising 
because a form of the inclusions of the second phase in this model has com- 
pletely different, layered, structure. In the low concentration expansion as well 
as directly in formula (64) one can see again that the divergencies appear in the 
limit z — *■ 1. 

(c) In case of almost equal phase concentrations x = 1/2 + e, e^Cl one 
obtains 

f(e,z)~l + 2ez, <r e// (e, z) ~ a + y/l - z 2 (1 + 2ez) . (68) 

The Keller - Dykhne formula (22) is reproduced for equal concentrations. 
The antisymmetric part of / has the following form 

- f 2 (e,z)-l 2ez 

fa{e, z) = —j- — = 1/2 - (69) 

2j(e,z) [1 - 4e 2 z 2 ] ' 

It follows from formula (41) that for e <C 1 and (or) z <C 1 the odd part f a 
coincides in the first order with the corresponding expression from the effective 
medium theory. The corresponding function $ is 



\/l — z 2 

In other words it is an analytic function of tz near ez — . Basing on this 
formula one can conjecture that an exact expression for / a (e, z) of this model 
will have a similar structure in general case with cr,; ^ or 1 — z ^ 1 . One 
must note that at the same time the formula (64) does not satisfy the equation 
(56) except of the trivial case x\ = X2- 



7. Comparison of different approximations 



For the comparison of the different expressions of the effective conductivity 
we have constructed 3D plots of f(e, z) in the EM approximation, in the FMSA 
approximation and for the random "parquet" model in FMSA- like approxima- 
tion (fig. 4, 5a). For a clarity they are represented at z < 0(cti < 02). 

It follows from fig. 4 and fig. 5a that all three formulas for a e , despite of 
their various functional forms, differ from each other weakly for z < 0, 8 due to 
very restrictive boundary conditions (25') and the exact Keller-Dykhne value. 
This range of z corresponds approximately to the ratio o^/ci ~ 10 _1 . For the 
smaller ratios a difference between these functions become distinguishable. 

It is useful also to compare a e in the FMSA approximation with cr™ of 
the Rayleigh model (RM), constructed numerically 14 j. Fig. 5b shows / in the 
FMSA approximation as a function of x and h = Trf for the direct comparison 

____ L-\-Z 

with the numerical results for RM from [14]. a r e m satisfies the self- reciprocity 
relation (17), but does not satisfy the duality relation (24). For this reason 
cr™ at adjoint concentrations are not related. The critical concentration x c 
of the RM in the percolation limit z — * 1 is x c = l- 7r/4w0.215 [T3|. It 
differs strongly from the critical concentration of self-dual systems x c = 1/2. By 



21 



compairing fig. 3 from ^3] an d our figs. 4b, 5b one can see that for z < 0.8 both 
<j e are very similar and close to each other. For z > 0.8 they begin to differ, our 
Ce becomes smaller than a^" 1 for concentrations in the range 0.2 < x < 0.8 due 
to a fact, that a e from (58) must follow the exact value (22) at x = 1/2. At the 
same time both conductivities remain close for concentrations in small regions 
near x = 1 and x = 0. This behaviuor can be explained by the difference in 
the geometry and symmetry of both systems, which becomes significant in this 
range of concentrations. For this reason one can expect, that for the RM with 
round inclusions deformed in a such way that the approximate symmetry (20) 
appears, its a r e m will become more similar to a e from (58). 

8. Percolation limit 

Now let us consider in the more details the derived formulas for <r e (e, z) 
in case when gi — ► 0(z — ► 1). It is clear that for regularly inhomogcneous 
medium one can always construct such distribution of the conducting phase 
that <7 e (e, 1) will differ from zero for all 1/2 > e > —1/2 . But in the case of 
randomly inhomogeneous medium the limit er 2 — > is equivalent to the well 
known percolation problem [17]. In terms of z it corresponds to the limit z — > 1 
and is also similar to the ideal conductivity limit a\ — > oo . Strictly speaking, 
an implementation of the duality transformation (14) is not obvious in this case. 
However, if one supposes that the dual symmetry relation (24) fulfills in this 
limit too due to a continuity then it follows from (24) that 

a e (e)a e (-e) = 0. (71) 

The relation (71) does not contradict to the known basic results of the perco- 
lation theory that <r e (e) = for e < and er e (e) ^ for e > 0. More- 
over it follows from the general formula (31) that in this case a s — \a a \ and 
a e (e) = cr a + \a a \. It gives 

i \ f 0, e < 0, 
^ = ( 2* a , e > 0. ^ 

This means that a behaviour of cr e (e) in the percolation theory is completely 
determined by its odd part. From a general discussion of the functional structure 
of the conductivity of the two-dimensional randomly inhomogeneous systems 
and experimental and numerical results it is known that in the percolation limit 
the effective conductivity a e must have a nonanalytical behaviour near the 
percolation edge x c = 1/2 or at small e > ^3 

<r e (e) ~ a x (x - x c Y ~ (Tie*, (73) 

where a critical exponent of the conductivity t is slightly above 1 and can be 
represented in the form t = 1 + 5. Since the values of this exponent found by 
the numerical calculations are confined to be in the interval (1,10 - 1,4) [17], 
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then 6 have to be small and belongs to the interval (0.1 - 0.4). It follows from 
general formula (72) that one must have 



z->l,e-(-0+ 



(74) 



It means that the function $(e, z) from (30) has at small e some crossover on 
z under z — * 1 from a regular (analytical) behaviour to a singular, powcrlikc 
one. At the moment an exact form of this crossover is unknown. For example, 
it can be of the form 



where $o(z 2 ) — ► 0, <&\{z 2 ) — » $i ^ and ^6^0 when z — > 1. Using 
(75) together with (30) and (32) one can show that this crossover takes place at 
1 — z ~ e 2t in agreement with result obtained from the similarity hypotheses 

on mi- 
it follows from the formulas for er e obtained in the sections 4-6, that one 
gets always a e — > in the limit (T 2 , except the region near x = 1. It 
means that these formulas obtained in the FMSA approximation are not valid 
in the limit 02 — > 0. This is confirmed by the appearance of the divergencies 
in the expansions of cr e in small concentrations in the limit z — > 1. This 
fact is a consequence of the assumptions made in the used approximations. 
For example, in case of the model of the layered (or striped) squares, where 
a deviation from the percolation type behaviour in the limit z — > 1 is more 
obvious, this is due to the "closing" (or "locking") effect of the layered structure 
in the adopted approximation in the limit oi — * 0. In order to obtain a finite 
conductivity in this model above threshold concentration x c one needs to take 
into account the correlations between adjacent squares. It is easy to show that 
near the threshold an effective conductivity is determined by random conducting 
clusters formed out of the crossing random layers from neighboring elementary 
squares. As is well known, the mean size of these clusters diverges near the 
percolation threshold [U| and for this reason the FMSA type approximation 
cannot be applicable for the description of er e in the limit z — > 1 (e > 0) 
and of the percolation problem. The divergencies of the expressions in the 
percolation limit mean strong changes of these expressions, which must result in 
a formation of the (widely believed) universal critical behaviour (73). It follows 
from figs. 4, 5 that the EM approximation overestimates a e , and both other 
formulas underestimate it in the region z — > 1, e > 0. We hope to investigate 
this limit in detail in other papers. 



Thus we have studied the possible functional forms of the effective con- 
ductivity of the two-phase SD system at arbitrary values of concentrations. A 
new approximate functional equation for a e , applicable for inhomogeneous self- 
dual systems and generalizing the duality relation, was deduced in the FMSA 




(75) 
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approximation and its solution was found. We have constructed also a hierar- 
chical (" parquet" ) model of the random inhomogeneous stripe type medium and 
have found its effective conductivity in the similar approximation at arbitrary 
phase concentrations. All expressions for the effective conductivity have the 
different functional forms. They 

(1) satisfy the self-dual symmetry and all necessary inequalities, 

(2) reproduce the general formulas for o~ e in the weakly inhomogeneous case. 

All these results confirm a conjecture that, in general, a e already of two- 
phase self-dual systems may be a nonuniversal function and can depend on some 
details of the structure of the randomly inhomogeneous regions ■ Analogous 
conclusions for SD systems with N = 3 phases were done earlier in the paper 
|22| . where a possibility to find a generalization of the exact Keller-Dykhne 
formula for case N = 3 was investigated numerically. 

The obtained explicit formulas can be used for an approximate description 
of the effective conductivity of some real (as regular as random) inhomogeneous 
systems like a corresponding formula of the EM approximation. Another im- 
portant property of these formulas is connected with a fact that they admit a 
generalization on a case of heterophase systems with arbitrary number of phases 
N (see [23 and my forthcoming paper). At the same time the formulas (58), (64) 
conserve their simple form for arbitrary N , while the EMA expression for o~ e 
becomes very complicated function even for N = 3. 
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(a) (b) 
Figure 4: Plots of f(e, z) in : a) the EM approximation, b) the FMSA approximation. 
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